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Abstract 

We apply a new coordinate space method for the evaluation of lattice Feynman 
diagrams suggested by Liischer and Weisz to field theories in two dimensions. 
Our work is to be presented for the theories with massless propagators. The 
main idea is to deal with the integrals in position space by making use of 
the recursion relation for the free propagator G{x) which allows to compute 
the propagator recursively by its values around origin. It turns out that 
the method is very efficient and gives very precise results. We illustrate the 
technique by evaluating a number of two- and three-loop diagrams explicitly. 


1 Introduction 


It is generally believed that QCD describes the strong interaction of quarks and 
gluons. QCD is, however, very complicated to investigate. Therefore, it is useful to 
study at hrst similar, but simpler models and apply the related knowledge gained to 
the real situation. The most frequently studied toy models for this 3+1 dimensional 
theory are the 1+1 dimensional models, e.g. the non-linear 0(n) a model in two 
dimensions, which have also found wide applications in the framework of the lattice 
held theory. They have not only the simple theoretical structure, but also are, 
due to low dimensionality, simple to simulate numerically. Another high interest 
in these two-dimensional models is to be found in perturbation theory which plays 
an important conceptual and practical role in lattice held theory. However, the 
evaluation of lattice Feynman diagrams is, in general, not easy. On the lattice the 
integrands become non-trivial function of internal and external momenta, which 
makes it difficult to apply the standard tools in continuum perturbation theory. 

A few years ago, Liischer and Weisz suggested a new position space method for 
the evaluation of lattice Feynman diagrams in four dimensional Yang-Mills theories, 
inspired by the observation of Vohwinkel [Q that the free lattice propagator can 
be calculated recursively by its values around origin. Their method concerns a 
technique of computing Feynman integrals in the coordinate space which allows to 
obtain very accurate results with little ehort and small amount of computer time. 
In our work on the computation of four-loop (3 function in the two-dimensional non¬ 
linear 0(n) a model [^, we have made extensive use of this technique. Since we have 
found it very useful, we would like, in the present paper, to describe the method in 
the case of two dimensions. This will be done by choosing some typical diagrams 
and evaluating them explicitly. 

To illustrate the basic ideas to be discussed in the present paper, we take the 
following 3-loop integral into consideration which will be treated in detail in the 
main text:Q 



t = —q — r — s, (1.1) 


where 



( 1 . 2 ) 


and p = {po,pi) denotes the line momentum. 

There are a few known conventional computer programs available to compute 
this integral. The numerical accuracy which can be achieved by these methods is, 
however, strongly limited and thus they are not very useful if one requires to have 
results with high precision. Considering the integral in position space, however, 
provides us with a new possibility to compute it very accurately with high efficiency. 

^Throughout the whole paper, we will work on a two-dimensional square lattice A and keep the 
lattice constant a being equal to 1. 
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In position space the integral (|1 . 1|) has the form 




(1.3) 


xGA /i=0 


where denote the lattice derivatives and G{x) is the position space free massless 
propagator 





(1.4) 


The free propagator G{x) diverges logarithmically at large x. After differentia¬ 
tion, d^G{x) falls off like and the snm (|1.3|) is therefore absolutely convergent. 
In principle, one can then compute Bi by summing over all points of x up to some 
hnite value, e.g. |x| < 30. However, since the convergence is slow, one requires to 
sum over large range of x if one would like to have precise result. Taking the limited 
computing availability into account, this is, of course, not the most efficient way. 
By studying the large x behavior of the propagator G'(x), it is however possible to 
improve on the accuracy systematically until the desired level of precision is reached. 
It is the aim of the present paper to show how this can be done. 

In next section we discuss the properties of the free propagator G{x) and, in 
particular, derive a recursion relation which allows us to express the propagator as 
a linear combination of its values near origin. These results will then be applied 
to the computation of lattice sums in section ^ where we take various examples to 
illustrate the procedure. It turns out that the method is quite efficient and gives 
very accurate results. Finally, in section we present a little more subtle case of 
evaluating Feynman diagrams with non-zero external momenta, particularly in the 
continuum limit. The usefulness of the position space method in this case as well 
is demonstrated by dealing with a number of diagrams explicitly. We conclude the 
paper by including some detailed computations in a series of appendices to make it 
more readable. 

2 Investigation of the free propagator 
and zeta functions 

In this section we prepare for the main work in the next two sections where we 
show how the position space technique can be applied to the evaluation of Feynman 
diagrams. This requires the study of large x behavior of the free propagator. We 
also need to derive the recursion relation for it which is essential to the computation 
of the integrals. 
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2.1 Properties of the propagator G{x) 

The propagator G{x) is a Green function for A, i.e. it satisfies the laplace equation 

- AG(a;) (2.1) 

where is the lattice h-function and A denotes the lattice laplacian 


A = E d;d,. 

/i=0 


( 2 . 2 ) 


Here, and 9* are the forward and backward lattice derivatives respectively: 

df,f{x) = f{x + fl) - f{x), (2.3) 

d*J{x) = f{x) - f{x - fi). (2.4) 

We note the relation d* = —for the backward lattice derivative, where is the 
adjoint of the forward lattice derivative dehned with respect to the inner product 
= {f,d^g). 

A rigorous derivation of the large x behavior of G{x) shows that in the limit 
X —> cxo the propagator diverges logarithmically and has the form 

G'^(x) = —^(Inx^ + 27 + 31n2) (2.5) 

4:71 

which corresponds to the leading behavior in the continuum limit. By making use of 
the laplace equation (|2.1|) , one can also work out the subleading terms systematically 
to obtain the asymptotic expansion. Up to order |x|“®, we hnd 


G{x) = G"(x) 


where 


Qi — 

Q2 = 




371 

'TE' 


„2\4 


^ 4523 

Qs = —r—[X 


2\6 


47r 

L(x2)3 

4 

-X , 


, 47 4 . , 

H-x^(x 

5 

7657 

— - 

-X I 


Qi + 


X 


2'|6 


Q2 + 


X 


2',9 


Qs + ■ ■ ■ 


20 


„4\2 


3716, 


„ 4 ^ 2 ,',„ 2^2 


2240 


56 ' 21 " ^ ^ ^ y ' 9 

In the above expansions we have introduced the shorthand notation 


X 


4n3 


x- = J2(xG^ 


fj,=0 

We remark that in two dimensions the general relation 

1 


2,2n+2 ^ ^2^2n _ tx^n-2 


f 2\2 4 

(X ) — X 


n > 1 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 
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is valid so that the expansion terms can always be expressed by the powers and 
since higher powers are to be reduced to these two powers (note x^ = 2). For 
example, from Eq. ( |2.11|) we hnd the identities 

‘ ( 2 . 12 ) 


3x'^x‘^ 


2\4 


X = - 

2 L 

X® = ^ {x'^Y + 2(x"‘)"*x^ — (x^) 
x^° = -x‘ 0(x‘)‘ — (X" 


„ 2\4 


(2.13) 

(2.14) 


The detailed discussion of the asymptotic behavior of G(x) is presented in Ap¬ 
pendix 

Now we would like to derive a recursion formula for the free propagator G{x) 
which makes it possible to determine G{x) for large x by its initial values around 
origin. The key relation for obtaining the formula is the identity, hrst observed by 
Vohwinkel §, 

G(x + fi) — G{x — fi) = x^if(x), (2-15) 

where H{x) is independent of fi and has the form 

d^p 


H{x) = 


Inp . 


(2.16) 


J-, (2ir)2 

By summing over p and using Eq. ( | 2 . 1 | ) for x 7 ^ 0, one can eliminate G{x + fi) in 
Eq. (|2.15|) and express H{x) like 


H{x) = 


ELo X. 


2G{x) — ^ G{x — fi) 

fj.=0 


(2.17) 


If one inserts this formula into Eq. ( |2.15|) , one hnally obtains the recursion relation 

2x„ 


G(x + /i) — 


ELnX, 


^/i=0 

which is valid for x 7 ^ 0 . 

Since the propagator has the properties 


2G(x) — ^ G{x — fi) + G{x — fi) 

/i=0 


(2.18) 


G(xo,xi) = G(xi,xo), (2.19) 

G(xo,xi) = G(-xo,xi), ( 2 . 20 ) 

G(x) is completely characterized by xq > Xi > 0. In this sector, Eq. (|2.18|) is a 
recursion formula which allows us to express G{x) as a linear combination of three 
initial values 


G(0,0) = 0, 
G(1,0) = 

G(l,l) = 


1 

■ 4 ’ 

1 

5 

TT 


4 


( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 



























i.e. the values of the propagator at the corners of the unit square. Therefore, the 
propagator has the general form 


G{x) = -\ri{x) - -r 2 {x), 
4 71 


(2.24) 


where ri{x) and r 2 {x) are rational numbers which are to be computed recursively 
by Eq. ( P.18|) . 

In this place we would like to stress the fact that the analytic knowledge of the 


initial values ( 2.21 )-( 2.23 ) is crucial to the precise determination of the propagator 
G{x) for large x. The reason for this is that ri{x) and r 2 {x) increase exponentially 
with X while G{x) grows only logarithmically and thus there is a large cancellation 
of terms on the r.h.s. of Eq. (|2.24|) . 


2.2 Properties of the function G 2 {x) 

Another function which often appears in the procedure of evaluating diagrams is 




2^2 


Goix) = 


iP 


l-n (27r)2 

The laplacian of this function satishes 

AG 2 ix) = -G{x) 


2\2 


(2.25) 


(2.26) 


One can also compute the function G 2 (x) analytically in the limit of large x 
and expand it systematically by using the expansion ( ^.6|) and the laplace equation 
(p.26|) . Up to order we hnd in this way 


G2(x) = 


X Ri + T 


where 


IGtt . 


Ri = lna;^ + 27 + 31n2-2, 




R2 

Rs = 
Ra = 


-(27 + 31n2), 


1 2 1 
= —Ina: H— , 

2 3 (a;2)2 

19 , 2x4 19 4 / 2\2 I 4\2 

= (a: ) + ’ 


, 4\3 


1680 


420 


21 


9 


X ) 


(2.27) 

(2.28) 

(2.29) 

(2.30) 

(2.31) 


For the explicit derivation of the large x expansion we refer to Appendix p. 

Similarly to the propagator G(x), one can also derive a recursion relation for 
G 2 {x). The hrst step for the derivation is the observation 


6 *2(x + p) - G2 {x - fi) = -x^H2{x), 


(2.32) 
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where H 2 {x) is independent of n and has the form 


H2{x) = G{x) + 


dvr 


(2.33) 


Eliminating G 2 {x + fi) in Eq. ( p.32D by using Eq. (|2.26|) gives the formula 


H2{x) = 


G{x) — 4:G2 {x) + 2 ^(0 G 2 {x — fi) 

fj.=0 


which leads straightforwardly to the recursion relation 


X 


G2{x + /i) — - , 


for X 7 ^ 0. 

Like the propagator G{x), due to the properties 

(j'2(Xo,Xi) = G2{Xi, Xq), 

G2{xq,Xi) = G2{-Xo,Xi), 


(2.34) 


^ G{x) — 4:G2 {x) + 2 ^(0 G 2 {x — /i) + G 2 {x — fi) (2.35) 

fj,=0 


(2.36) 

(2.37) 


G 2 {x) is completely characterized by Xq > Xi > 0. In this sector, the recursion rela¬ 
tion (|2.35|) allows G 2 {x) to be expressed in terms of G{x) and a linear combination 
of its values at the corners of the unit square, 


^2(0,0), (^2(150), G2(1,1). 

These initial values can also be computed analytically with the results 


(2.38) 


^2(0,0) = 0, (2.39) 

^2(1,0) = 0, (2.40) 

0,(1.1) = T, (2,41) 

2.3 Zeta functions Z{s,h) 

In the procedure of evaluating the lattice sum like Eq. o. zeta functions appear. 
As a preparation for the computation of such lattice sums in next section, we discuss 
the evaluation of the zeta functions Z{s, h) associated with the square lattice A. 
Let us dehne a genealized zeta function through 

Z{s,h) = ^'h(x)(x^)“", (2.42) 

where h{x) is a harmonic homogeneous polynomial in the real variables x = (xq, Xi) 
with even degree d and s complex number with Res > 1 -|- d/2. The primed 
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summation symbol implies that the point a: = 0 should be omitted in the summation. 
Note that the sum ( 2.42|) is absolutely convergent in the specihed range of s. 

The zeta functions can be evaluated as follows. We hrst introduce the heat kernel 


k{t, h) = ^ h{x)e 


— TTtX^ 


(2.43) 


xGA 


In terms of the heat kernel, Eq. (p.42|) can be rewritten in the form 


Z(s, h) = 


IT 


dtt^ ^[k{t, h) — h(0)]. 


r(s) Jo 

Now we observe that the heat kernel has the property 

k{t, h) = {-ly/H-'^-^kil/t.h) 

which follows from the poisson summation formula 

—iqx^—Tvtx^ _ 1 ^-{q+2TTxf/{4TTt) 


e "'^"'e 

xeA 


xGA 


(2.44) 


(2.45) 


(2.46) 


and the condition that h{x) is a harmonic function. 

If we divide the integral (p2.44| ) into two parts and make use of the identity {\2,A5\ ), 
we hnally obtain the representation 


Z{s,h) = 


TT 


h{0) 


+ 


r(s) l s{s — 1) Ji 


dt 


\k{t, h) - h(0)l 1 . (2.47) 


In the integration range 1 < t < cxo, the large x terms in the series (|2.43|) are 
exponentially suppressed. Eq. (|2.47| ) can therefore be computed numerically very 
precisely with little effort, e.g. by using symbolic language MATHEMATICA. 


3 Evaluation of lattice sums 

By dealing with a number of integrals explicitly, we would now like to show how 
the results of the last section are to be applied to the evaluation of lattice Feynman 
diagrams. Concretely, we will consider in the following the two- and three-loop 
diagrams of hgure |^. 


3.1 Two-loop Feynman integrals 

As a hrst example let us consider the integral, encountered in the 3-loop (5 function 
computation in the 0{n) a model. 


Ai — 


/- d^k dH 1 

(27r)2 (27r)2 ^o ^ 2 ’ 


s = k + I 


(3.1) 
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Figure 1: Two- and three-loop diagrams 


which corresponds to the hrst diagram in hgure 
In position space this integral can be written as 

where 


(3.2) 


^i(^) = '£(9,- d;)G{x){[(d^ - d;)G{x)f + [(8„ + 0;)G(a:)]h (3.3) 

/^=0 

G{x) diverges logarithmically at large x and an explicit computation of large x 
behavior of Eq. (|3.3| ) shows that Ai{x) falls off like \x\~^. The sum ( p.2|) is therefore 
absolutely convergent and a hrst approximation to Ai may be obtained by summing 
over all points x with (say) |x| < 30. A direct summation in this way gives a result 
as a hrst approximation 

Ai = -0.04166(1). (3.4) 

We can now improve on the accuracy by investigating the asymptotic behavior 
of Ai{x). From the asymptotic form of G{x) [Eq. (|2.6|) ], we deduce that 


where 


{^i(^)}, 


27r^ 


X' 


2\2 


as 


= {ii(x)}_^ + 0(|x|-«), 

(3.5) 

x^ 10 x'^ ix^A ^ 

^(3.2)4 (^2)3 1 ^“(^2)5 

|. (3.6) 

of Eq. (|3.6|) into consideration, Ai can 

be rewritten 

= ^ + ^i(O) ) 

(3.7) 


(3.8) 

tcEA 


(3.9) 


’>=-? 

(3.10) 


where 
















The sum where the asymptotic series has been subtracted converges much 
more rapidly than the original sum Ai. This subtracted sum can therefore be com¬ 
puted very accurately by summing over all points x with |a:| <30. As for the term 
A^ from the asymptotic expansion, it can be computed with help of zeta functions 
very accurately as well. For the purpose of applying the zeta functions, we transform 


the asymptotic expansion ( p.6|) in the following way: 




2 1 
+ 


°° (27r)^ I 


ho{x) + 


2 4 

+ 


L(a:2)4 (a.2)5j 


hi{x) 7 y^h 2 (a;) 


[X 


(3.11) 

Here, h^^x), hi{x) and h 2 {x) are the harmonic homogeneous polynomials given by 


ho{x) = 1, (3.12) 

hi{x) = — 3(x^)^, (3.13) 

h2{x) = ‘i2{x‘^Y — A&{x‘^Yx‘^+ VJ{x^Y. (3.14) 


In Appendix y we discuss how one computes the harmonic polynomials. The sub¬ 
traction term A^ can then be written as a sum of zeta functions: 




( 2 ir) 


:Jy2Z (2, hg) -|- Z (3, Hq) -\- 2Z (4, hi) -\- AZ (5, hi) -\- “AZ (7, h2)| (3.15) 


to which now the solution (|2.47| ) can be applied. 

Summing all three contributions together, we get hnally 


yli = -0.041666666666(1). (3.16) 

It seems very likely that Ai is — Indeed, the analytic computation of Ai conhrms 
this hypothesis (see Appendix 0). We note that it is straightforward to determine 
Ai more precisely just by including more terms in the large x expansion of (|3.6|) . 

It is interesting to compare our determination with the earlier estimation by 
Weisz [H who states the value 


Ai = -0.0416666(2). 


(3.17) 


We clearly see that the result from our method is much more precise than that of 
Weisz. 

As a second example, we consider the following integral which also appears in 
the 3-loop {3 function computation (see hgure |l|) 


Ao — 


r d^k d^l s^-P-P P 


l-n (27r)2 (27r)2 (, 

Here, we made use of the notation of Eq. (|2.10|) . 


;2'i2 ■ 


S = k + 1. 


(3.18) 
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The evaluation of this integral proceeds in exactly the same way as that of Ai. 
In terms of the propagator G(x) and the function G 2 {x) [Eq. (|2.25|) ], in the position 
space the integral can be expressed by 

A2 = --^A2{x), (3.19) 

^ xga 


where 

A2{x) 



The terms in this sum fall off like |x| and a direct evaluation without subtraction 
gives 

A 2 = -0.18465(1). (3.21) 

The convergence can be improved by subtracting the asymptotic form of A 2 (x) 
which can be deduced from the asymptotic forms of G(x) and G 2 (x). Up to order 
we hnd for the large x expansion of A 2 {x) 


A2{x) = [A2{x)^^ + 0{\x\ ®), 


where 


1 

IGtt^ 


5 

(^ 


ho^x) + 



(3.22) 


(3.23) 


The subtracted sum 


^2 = - {^2(2:)}^ 


(3.24) 


is now rapidly convergent and may therefore be computed accurately. On the other 
hand, the subtraction term 

T = (3.25) 

xGA 

can be expressed as a sum of zeta functions: 

A 2 = (3, ho) + 8Z (4, hi) + 2Z (5, hi) + 25Z(7, h2)| (3.26) 


which can be computed very accurately, too. 
Since 


we obtain hnally 


A2 — —- ^42 + ^42+ ^42(0) 


A 2 = -0.1846545169(1). 


(3.27) 

(3.28) 


This result is also to be compared with the earlier estimation by Weisz [|] who states 


A 2 = -0.1846544(1). 


(3.29) 
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3.2 Three-loop Feynman integrals 

Now we come to the 3-loop integral of Eq. ( p..l|) which we have already discussed in 
Introduction (see hgure |l|). The summand of Bi from Eq. (|1.3|) 


B,{x) = -z [a,G( 

)j,=0 


X 


(3.30) 


falls off like |a:| ^ at large x and a direct summation without subtraction gives 


Bi = 0.0169611(1). 


(3.31) 


For the acceleration of the convergence, we expand Eq. (|3.30|) for large x. Up to 
order we hnd 

Bi(x) = {Bi(x)}^ + OHxr^), 


where 


From the identity 


X 


3 _ 6^1 + 16 '“^ 


4\2 


oo 327rn (a;2)4 (a;2)3 "(^2)5 "(a;2)7 

Bi = B[ + B'' + Bi(0) 


with 


xgA 

Bi( 0 ) = 

and the zeta function representation 
1 


b; = Z![Bi(x) - {b,(i)}_ 

xGA 

b'; = z'{b4x)} 


128 


(3.32) 

(3.33) 

(3.34) 

(3.35) 

(3.36) 

(3.37) 


B'l = 


647r^ 


|3.Z(2, ho) -|- 4.Z(3, ho) hi) -|- 9.Z(5, hi) Z(7, h2)|, (3.38) 


it is straightforward to compute Bi to get a more precise result 

Bi = 0.016961078576(1). 


(3.39) 


We remark that this integral has been computed by Caracciolo and Pelissetto 
who state the value 

Bi = 0.016961. (3.40) 

As another example, we consider the 3-loop integral of the form 


Bo = 


n cPq d?r (Ps ( D)l.=o 
U (27r)2 (27r)2 (27r)2 ^ 2 ^ 2 ^ 


t = —q — r — s. 


(3.41) 
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In terms of the position space propagator G(x), the integral can be expressed by 

(3.42) 


B2 = Z^2(x), 

xeA 


1 1 


(3.43) 


where 

^2(x) = J2 12 [d^dvG{x) 

^=0u=0 

The terms in this sum fall off like |x|“® and converge very well. Therefore, the direct 
summation already gives a quite precise result: 


B 2 = 0.13661977236(1). 


(3.44) 


Through exactly the same procedure as Bi, we obtain the following identity for 
B 2 which should give more accurate result: 

B 2 = E'[^ 2 (x) - {52(x)} J + ^{ 3 ^( 4 , ho) + Z(8, h 2 )} + 5^(0) (3.45) 


cGA 


with 






87r‘^l(x^)^ (a;^)® (a;^)®]’ 

~ , 34 36 15 3 1 

^ 2 ( 0 ) = 

TT^ TT-^ TT^ TT 4 


(3.46) 

(3.47) 


The numerical evaluation of Eq. ( |3.45|) then yields a more precise value 

B 2 = 0.136619772367581(1). (3.48) 

This number can also be compared with that by Caracciolo and Pelissetto P: 


Bo = 0.1366198. 


Finally, we consider the 3-loop integral of the form 


(3.49) 


S, = 


d^g dV d^s {q + r) - f {s + t) - P 


l-n (27r)2 (27r)2 (27r)2 

t = —q — r — s. 


The position space representation of this integral is given by 

1 


where 


S3=T i: B,(x), 

^ xeA 


Bsix) = { E [((^M - 9*^)G{x)y + + d*^)G{x)) 

^ 11=0 


2i 


(3.50) 


(3.51) 


(3.52) 
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The evaluation of the sum (|3.51| ) runs in exactly the same way as above and we 
get the result 

S 3 = 0.09588764425(1), (3.53) 

where we have evaluated the sum by subtracting the asymptotic terms from the 
summand in Eq. (|3.51|). Finally, this result is to be compared with that in Ref. 1^: 


S 3 = 0.095887. 


(3.54) 


4 Evaluation of diagrams with external momenta 

In this section, we consider the Feynman diagrams in hgure ^ with external mo¬ 
menta p. We would like to evaluate the integrals in the continuum limit which 
corresponds to taking p to zero since we are using lattice units and there are no 
mass parameters. In other words, we have to work out the leading terms of the 
integrals in an asymptotic expansion for p — 0 . 



Figure 2 : One-, two- and three-loop diagrams with non-zero external momenta 

For the evaluation of diagrams with external momenta p, it is necessary to solve 
the p-dependent zeta functions. As a preparation for the actual treatment of such 
diagrams, we therefore discuss at hrst the zeta functions which are depending on 
the external momenta. 

4.1 Zeta functions s) 

Let us consider the p-dependent zeta functions of the kind (s > 1) 

„-ipx 

= (4-1) 

With help of the heat kernel [Eq. ( |2.43|) ], we can rewrite Eq. ( |4.1j) in the form 

z(p, s) = / dt f^ . ( 4 . 2 ) 
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We are interested in the evaluation of this integral in the limit p —0. For this 
purpose, we consider the integral by dividing the integration range into two parts 


Z{p,s) = Z'{p,s) + Z"{p,s), (4.3) 

where 

Z\p,s) = -l), (4.4) 

Z"{p,s) = ^ -l). (4.5) 

r(s)./o ^ 


Now, note that Z'{p, s) is well-dehned for all s > 1. We can therefore set p to zero 
in that function. Then, Z'(0, s) can be evaluated numerically very accurately since 
the large x terms in the series of the function are exponentially suppressed. 

On the other hand, Z"{p, s) contains singular term for small p. For the evaluation 
of this function, we hrst apply the Poisson summation formula [Eq. (|2.46|) ] to rewrite 


Z^P.s) = 


TT 


r(s) Jo 


r dt t‘-' 

^-1 ^ g-(p+27ra;)2/(47rp _ ^ 

Jo 

xGA 


(4.6) 


In the series of the integrals appearing in Eq. ([4.6|) , the only dangerous term for 
p —0 is that with x = 0. Hence, we split the integral Z"{p, s) into the regular and 
singular parts 

Z"{p, s) = Z"{p, s) + Z''{p, s), (4.7) 

where 

z;'(p,s) = /°°dfr*e-(^’+2-0W(40^ (4.8) 

2:{p,s) = ^ (4.9) 

1 (s) Jo 

After having set p to zero in the well-dehned function y (p, s), the integral Z"{0, s) 
can also be calculated numerically with high precision. As for Z"{p, s), this function 
can easily be evaluated analytically in the continuum limit by appropriate separation 
of the singularity. 

Summarizing the whole terms, we get hnally for p ^ 0 



with p = Att/p^. 
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4.2 Scalar one-loop integral 

We consider the integral 


DM = / 


(ik 


-n (27r)2 


s = k + p 


(4.11) 


which corresponds to the first diagram in fignre ^ By applying the coordinate space 
method, we wonld like to show that -Di(p) has the form 


-Di(p) = Cl Inp^ + C2 + O(p^) 

and determine the coefficients ci and C 2 .[| 

In position space, Eq. ( 4.11|) has the form 

For large x, the fnnction AG(x)^ has the asymptotic behavior 

AGixf = {AG(a:)2}^ + 0{\x\-^) 

with 

{AG(x)q 


1 1 


4 7 1 X 


' oo 27r2 x"^ 

This snggests the splitting of Di{p) into fonr parts: 


2x‘^ 


where 


DM) = -\d[{p) + DM) + D';'(p) + AG(x)'y.„ 


DM) = - {AG(i)q^ 


cGA 


DM) ^ 


1 1 


cGA 

D';'(P) = 

xGA 

AG(a;)^|a;=o = t- 


27r^ x‘^ 


/e 


—tpx 




(4.12) 


(4.13) 

(4.14) 

(4.15) 

(4.16) 

(4.17) 

(4.18) 

(4.19) 

(4.20) 


The fnnctions D[{p) and 77'/(p) are well-defined and therefore the limit p —0 
exists. After having set p to zero, 71((0) can be compnted nnmerically with high 
precision by applying the position space method. We receive the resnlt 


L)((0) = -0.14215315037(1). 


(4.21) 


^In this section 0(p") stands for a remainder R{p) such that limp_,o i?(p)/|_p|" '^ = 0 for all 
e > 0. 
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D"( 0 ) can also be computed very accurately by making use of the solution of zeta 
functions [Eq. ( p.47 )]. A straightforward calculation gives 


L)"(0) = 0.275883297410294(1). 


(4.22) 


The function Zl("(p), however, diverges in the continuum limit, which makes it neces¬ 
sary to apply Eq. ( [4.10| ). After some algebra and numerical evaluations, we hnd 


D'l'ip) = -lnp2 + 0.1678588533415551(1). 

27r 


(4.23) 


Summing all four contributions together, we conhrm that Di{p) has indeed the 
form of Eq. ([4.12|) , and for the constants Ci and C 2 we obtain 

1 


Cl = 


27r’ 


C 2 = -0.55158900038(1). 


(4.24) 

(4.25) 


4.3 Scalar two-loop integral 

Now let us consider an example of two-loop diagrams. Concretely, we would like to 
discuss the following integral shown in the second diagram of hgure ^ 


D2{p) = 


r- dH d^s {k + l) -P- f 1 

L (^ (^ (^ pp 


^<^p-i-fc-i-z+s,o *^p-i-fc-i-z,o) • (4.26) 


We have used here the shorthand notation Spp = (27r)^5*^^)(p). 

In position space, the integral (|4.26|) can be written in the form 

D2(p) = -) E 

^ xgA 


with 


D 2 W = gwj:{[(5„-5;)g(i)] +[(5„ + s;)gw] }. 

n=o 


(4.27) 


(4.28) 


For the evaluation of the sum ( [4.27] ), we hrst observe the expansion terms of D 2 {x) 
for large x. Up to order we hnd 


where 




D2(x) = {r>2(l)}_^ + 0(|l|-®), 
1 r In 27 -|- 3 In 2 


(4.29) 


47r^ I 
7 1 


x^ 


+ 


x^ 


2 (a;2)2 L 


+5- 


X 


(X 


2'|4 


In -I- 3 In 2 -I- 27 - 

7. 

9 2 

In -|- 27 -|- 3 In 2- 

15 


(4.30) 
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Now we separate the singular part in Eq. (^4.27|) and rewrite it as 


where 


D2(p) =+ D'^ip) + d:^'{p) 


DW = J2'o“‘"[D2(x)-{D2(x)} 

x£A 

d'Hp) = 


cSA 


dvr^ V 


+ 






D^P) = 


1 


dTfS 


x&A 


ipx 


In x'^ 27 + 3 In 2 


+ 


x^ 


x^ 


(4.31) 

(4.32) 

(4.33) 

(4.34) 


We note that £> 2 ( 0 ) = 0. 

The function Di^lp) is well-dehned in the limit p ^ 0 and can thus be computed 
numerically in this limit by applying the position space method. We obtain 


L)^(O) = 0.0705019591(1). 


(4.35) 


112 ( 0 ) can also be evaluated very accurately by using the solution of zeta functions 
[Eq. (1^) ]: 

D”{0) = -0.1385014044234234(1). (4.36) 

The function D 2 {p), however, diverges in the continuum limit. To this divergent 
function we can again apply the solution ( |4.10|) as in the case of the one-loop integral. 
The only difference here is the appearence of a logarithmic term in zeta function 
which can, however, be evaluated by differentiating the appropriate zeta function 
with respect to s. In this way we obtain 


nW/' ^ ^ /'i 2^2 I 5 In2 

^2 [P) = ) + 


lnp2 _ 0.0841257673993860(1). 


87r2^ (27r)2 

Summing all three terms together, we hnally get for D 2 {p) 


D2{p) = 


-(Inp- 


2\2 


51n2 


with Ki given by 


(47r)2^ 87r2 

= 0.0760626064(1). 


Inp^ -|- Ki 


(4.37) 

(4.38) 

(4.39) 


4.4 Scalar three-loop integral 

As the last example, we consider the 3-loop integral of the form (see hgure 

d^k dH d^s dH 


Ds{p) = J_ 


-pr (27r)2 (27r)2 (27r)2 (27r)2 kH^H^ 

^p+fc+i+SjO ^p+k+l+t,0 T ^p+k+l,0 


(4.40) 
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In position space this integral can be expressed by 


^ xsA 


(4.41) 


where 


D,{x) = Gixf { [(Sm - 8;)G{x)] + [(8„ + d;)G{x)] }. (4.42) 

11=0 


Up to order |a;| we obtain for the large x expansion of the fnnction D 3 {x) 


D,{x) = {d,{x)}^ + 0{\x\-^), 


(4.43) 


where 




1 1 
1 


(Inx ) + 2Qo \nx + Qq 


x^ L 2 


^(Inx^)^ + (1 - 7Qo) Inx^ - JQq + Qo 


+ 


X 


[X 


2'|3 


5(lnx^)^ + (lOQo - g) Inx^ + 5Ql - -Qo 


(4.44) 


with Qo = 27 + 3 In 2. 

The snm ( 4.41 ) can now be evaluated in exactly the same way as D 2 {p). The 
only difference in this case is that, due to the quadratically logarithmic term in 
the expansion ( [4.44| ), it is necessary to differentiate zeta functions two times with 
respect to s. After a straightforward computation, we get for D^ip) 


Dsip) = 


with K 2 given by 


327r3 


-(Inp^)^ — 51n2(lnp^)^ + (5 In 2)^ In 
3 


K2 = -0.01802345(1). 


+ K2 


(4.45) 

(4.46) 
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A Large x expansion of G{x 


In this appendix we derive for the free propagator G{x) given in Eq. o the large x 
expansion of Eqs. (pA|) and (p.6|) . 
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As a first step, we use Cauchy’s integral formula to do the po integration and for 
simplicity set x = (a:o, 0 ). ThenQ 


fTT dp gXo In i [p2+2-y^p2(p2+4)] _ ^ 

G{xo) = -/ =-• 

Jo TT /p2(^p2 _|_ 4^ 


(A.l) 

Change variable y = sinp/2, then p = 2y, dp = 2dy /and Eq. (|A.1|) becomes 

1 /•! dv — 1 

g(^o) = ^ 4 > (A-2) 

ZTT JO ?/ V 1 ~ 2/ 

where fiy) is dehned by 


f{y) = In (l + 2 y^ - 2 yyjl + y^'j. 

Now, divide the integral of Eq. (|A.2|) into two parts 


G(xo) — — G(xo) + G(xo) 
2 tt l 


where 


G{xo) — 
G{xo) = 


ri dy r 


/o y 
/o y 




1 - 


VI - 1 /^ 


(A.3) 

(A.4) 

(A.5) 

(A. 6 ) 


We are interested in the evaluation of G(xo) in the limit Xq —>■ oo. For this 
purpose, we perform another variable change t = x^y in G{xq). In the continuum 
limit (xo —> oo), Eq. (|A.5|) then reduces to 


G{xo) 


r^o , e 
dt ■ 


-2t 


(A.7) 


The integrals of Eqs. (|A. 6 |) and (A-7|) can now be computed analytically and we 
hnally get for the leading terms in the continuum limit 


G‘’(xo) = - —( 21 nxo + 27 + 31n2). 
dvr 


(A. 8 ) 


Using rotation invariance in the continuum limit, the leading behavior for G{x) is 
then equal to 

G^(x) = —^(Inx^ + 27 + 3 ln 2 ). (A.9) 

4:71 


^For a given function F{x), the simplified notation F(xq) stands for the exact F{x)\x=(xo,o)- 
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The subleading terms can also be computed systematically if we make use of the 
laplace equation (El- We derive, at hrst, the hrst subleading terms through the 
ansatz 


G'(x) = + G'W(a^) + 0{\x\-^) (A.IO) 


with 


G-q.) = -- 


1 X‘ 

Cl —+ C2 


X^ 


2;2)3_ 


(A.ll) 


For determination of ci and C 2 , we evaluate the lattice laplacian of G{x) in Eq. (|A-10D 
to get 


AG'(a:) = 


6 


X 


(X 


2'|4 


+ Ci- 


[X 


2\2 


+ 12C2 


X 


4 n 


[X 


2\2 


X 


2'|4 


} + 0(|i|-'). (A.12) 


dvr I (x^)^ 

From the condition that the r.h.s. of Eq. (|A.12|) should be zero, we obtain 

1 2 


Cl = 


C 2 = 


(A.13) 


Now, we compute the next subleading terms and try with the ansatz 


G(x) = G%x) + + G^^\x) + 0(|x|-®), 


where 


G^^\x) = - 


dvr 


C3- 


[X 


2\2 


+ C4- 


X 




+ C5- 


[X 


A \2 


X2)6J 


The lattice laplacian of this function is given by 




d80- 


X 


[X 


2'|5 


- 320- 


(X 


4\2 


[X 


217 


+ C3- 


+C4- 


12 


^2)3 


C 5 


16 


[X 


213 


™4 

72™ + d8‘ 


(X 


215 


[X 


217 


16 

x2)3 

+ 0{\x\ 


Since the r.h.s. has to vanish, we obtain 

2dc3 + I 8 C 4 — 95 = 0 , 


C 5 = 


20 

y 


(A.14) 

(A.15) 


(A.16) 


(A.17) 


Unfortunately, the condition from the laplace equation alone is not enough to 
determine the coefficients completely and we need another constraint for them. It 
can be obtained by computing the subleading terms of G(xo) in Eq. ( |A.1|) to this 
order. After some algebra, we hnd 




1 1 




„-6l 


(A.I 8 ) 
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(A.19) 


By comparing this with Eq. ( |A.14|) for x = (xq, 0), we get the relation 

43 


C 3 + C 4 + C 5 — 


120 


From Eqs. (|A.17|) and (|A.19D, the coefficients C 3 and C 4 can now be determined: 


C 3 = 


371 


C 4 = 


47 


(A.20) 


The computation of the third subleading terms proceeds completely analogously. 
After a lengthy, but straightforward computation, we hnd 


G(x) = G\x) + + G^^\x) + 0(|xr®) 


with 


G^^\x') = - 


47r 


4523 1 


7657 x^ ^ 3716 (x^ 2240 (x^) 


4\3n 


56 (x^)^ 21 (x^)^ 7 {x‘^y 9 (x^)^ 


(A.21) 


(A.22) 


B Large x expansion of G 2 {x) 


We derive here the large x expansion of Eq. (|2.27|) for the function G 2 {x) given in 
Eq. (^) . 

We hrst do the po integration and for simplicity set x = (xq, 0). Then 


^ 2 ( 2 ^ 0 ) — 


dp r (xo + l)e*'’*“ 2 (p ^+2 _ 1 _ 


/O TT I 


(p2 + 2 - P)[e^oln2(P+2-P) _ 1] _ ^(^2 _ py 

P3 


(B.l) 


with P = ^Jp^ijP' + 4). By changing variable y = sin p/2, Eq. (p.l|) transforms to 

1 j-^ dp e-°/(")(2p^ + 2xopyr+^+ 1) - [1 + 2(1 - xg)p^] 2 ) 


Ibvr Jo 


(1 + 2/2)71^ 


with /(p) given in Eq. (A-3|) . 

Now, divide the integral of Eq. ( |B.2|) into two parts: 


^2(2:0) — T7P G2{xo) + G2{xo) 
Ibvr L 


where [F = (1 + y^y/T^^y^ 


G2(xo) — 


ni I QXofiy) (2^/2 P 2xoy^/TTf + 1) _ 




lo y^ { Y 

l + 2p2' 


(-) ^ (-f) 


(B.3) 


l + (l-2x2y ^,(b.4) 
1 


(B.5) 
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We are interested in the evaluation of G 2 {xo) in the continuum limit a:o oo. 
For this purpose, we perform the variable change t = xoy. In the continuum limit, 
Eq. ( [B.4|) then reduces to 

< 52 ( 3 ^ 0 ) ~ + + 

1 r^o 

+ - / dte“^‘(l + 2t). (B.6) 

O JO 

In the limit of inhnitely large Xq, these integrals can be evaluated analytically to 
yield 

< 52 (^ 0 ) ~ xl (21na:o + 21n2 + 27 - 1) - Inxo - ln2 - 7 + ^ • (B.7) 

o 

As for G 2 {xo), we can also compute the function analytically: 

G' 2 (xo) = Xq (In 2 - 1) - ^ In 2 . (B. 8 ) 


Summing two contributions together, we have finally for the leading terms in the 
continuum limit 


G *2 (^ 0 ) = 77 — xl (2 In xo + 3 In 2 + 27 — 2) 


Ibvr L 


In xn-In 2 

2 


In 

7 H— . 


(B.9) 


Unfortunately, from Eq. (P3.9| ) alone one can not derive the corresponding func¬ 
tion G^ix) unambiguously since G 2 {x) in the continuum limit is not rotation invari¬ 
ant. We therefore apply the laplace equation ( p.26| ) from which we obtain for G 2 {x) 
the following leading terms in the inhnitely large x limit: 




Ibvr L 


1 x^ 1 

x^ Inx^ -I- (3 In 2 -|- 27 — 2)x^ -|- In x^ + K 2 ■ (B. 10) 


3 (x2) 


Comparison of this result for x = (xq, 0) with Eq. (|B.9|) gives hnally 


= -^ , K2 = -^(31112 + 27) . 


(B.ll) 


The subleading terms for G 2 {x) can also be worked out systematically by apply¬ 
ing the laplace equation (p.26|) and the asymptotic expansion for G 2 {xq) of Eq. ([B.l|) . 
We do not present here the detailed computations leading to the result of Eq. ( p.27| ) 
since they are too long and also proceed very similarly to the case of G(x). 


C Derivation of harmonic functions 

We discuss here how one computes the harmonic homogeneous polynomials appear¬ 
ing in the main text. 
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Let us consider the following polynomial functions of degree 4n 


hn{x) 


n 

r=0 


(C.l) 


with y 0. We would like to derive the conditions for the coefficients so 
that hn{x) are harmonic functions. 

For determination of the conditions, we hrst note the relation 


A, 




= 4 


2r(r-l)(a:^)^-2(a:2)^+3 


+3r(2r — + s(4r + s){x‘^Y{xy 


4\r /^2\s—l 


(C.2) 


where Ac denotes the continuum laplacian Ylli=o 9“^/dx^ In deriving the relation 
above we made use of the identity (|2.12|) . By using Eq. (|C.2|) we obtain 


IL 

Achn{x) = 4y^ — Y) + + l)(2r + l)0{n — 1 — r) 


r=0 


2 A^+ 2 (^ + l)(r + 2)6{n — 2 — r) 


A\r { 2\2n—2r—l 
Jb 


(C.3) 


where 6{t) is a step function dehned by 


e{t) 


1 for t > 0 , 
0 for f < 0 . 


For hn{x) to be harmonic, the r.h.s. of Eq. (|C.3| ) should be zero, which leads 
to the conditions for the coefficients CY'* to satisfy. Now, we investigate these 
conditions by dividing the n + 1 terms in the summation contributing to A^hnix) 
into the following three cases: 

1. For r = n\ the contribution to A^hnix) is identically zero and there is therefore 
no constraint on the coefficients 

2. For r = n — 1: the contribution to Achn{x) is zero if 

3nCy^ + = 0 . (C.4) 


3. For r < n — 2: the contribution to Achn{x) is zero if 

4Cy\n^ - rY + 3C'i?i(r + l)(2r + 1) - 2Cl%{r + l)(r + 2) = 0 . (C.5) 


It is now straightforward to derive the harmonic homogeneous polynomials in¬ 
cluding those of Eqs. ( |3.12|) - (|3.14|) by solving the recursion relations ( |C.4| ) and (|C.5|) 
up to an arbitrary overall normalization factor. 
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D Analytic computation of Ai 


This appendix where we evaluate the two-loop integral Ai given in Eq. (O) ana¬ 
lytically is due to communication with A. Pelissetto. 

We first use the identity cos ^ = sin and the permutation symmetry bet¬ 
ween k, I and (—s) to transform the given integral to 


A — t; 


r (fk dH 


3J-n (27r)2(27r)2A;2/252^ 


Eh' 




.sms,, 


sin k^ — sin 


(D.l) 


11=0 


Now, we use the following trigonometric identity being valid for a -|- /3 -f 7 = 0 


sin a -|- sin /3 -f sin 7 


—4 sin 


tt . /d . 7 

— sm — sm — 
2 2 2 


(D.2) 


to rewrite 


Ai 


4 

3 



d^k dH 1 

(27r)2 (27r)2 


El 1 ^ 

y^s^sm —sm-sm—. 

/.t =0 


(D.3) 


Noting the relation sin ^ we see that this integral is, apart from a factor, 

same as A® in Appendix A of Ref. [^] which is known analytically. From this, we 
obtain 



(D.4) 
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